Comment on "Quantum back-reaction through the Bohmian particle 



L.L. Salcedo^ 

Departamento de Fisica Moderna, Universidad de Granada, E-18071 Granada, Spain 

(Dated: February 1, 2008) 



In a recent Letter [l| Prezhdo and Brooksby have 
forwarded a solution to the long-standing problem of 
properly dealing with a mixed quantum-classical system, 
namely, to couple the classical sector, not directly to the 
quantum one, but to its associated Bohmian particle 0. 
The precise proposal is to make replicas of the quantum- 
classical system, each one with a different position of 
the Bohmian particle obtained by sampling the "initial" 
quantum sector density. Each copy of the system is then 
to be evolved as follows: i) the quantum sector evolves as 
usual, i.e., by taking the classical degrees of freedom as 
time-dependent parameters of the quantum Hamiltonian, 
ii) the Bohmian particle evolves a la de Broglie-Bohm, its 
velocity depending solely on the quantum degrees of free- 
dom and iii) the classical sector evolves following Hamil- 
ton's equations but this time with the Bohmian particle 
position as time-dependent parameter in the Hamilto- 
nian. The expectation values of observables are obtained 
by averaging over the "final" sample of quantum-classical 
degrees of freedom. In practice this proposal coincides 
with that in |^ (in its simplest version) although there it 
was derived as an approximation rather than introduced 
as a postulate. In this Comment I want to point out some 
limitations of the proposed prescription if it is pushed too 
far. 

Without loss of generality, it can be assumed that 
the initial state is given by some probability distribution 
n{u, ip) defined on the combined manifold with coordi- 
nates u = {x, k) (the position-momentum of the classical 
particle) and (the normalized wave-function of the 
quantum particle). The full distribution f{u,'ip,y) in- 
cluding the Bohmian particle (with coordinates y) will 
then equal n{u, ip) |V'(y)P- In this way / is initially fully 
determined by n. The evolution equations then define a 
flow for / in the manifold {u,ip,y). The flow of / has 
no diffusion (an initial Dirac delta distribution evolves 
as a delta all the time). Of course, its marginal den- 
sity n{u,ip;t) ~ J dyf{u,^p,y;t) does show diffusion, 
which is interpreted as branching and back-reaction of 
the quantum-classical system. 

The first problem of the proposal is that the initial 
relation noted between / and n is not preserved by the 
evolution, that is, at later times the Bohmian particle 
is no longer distributed as |V'(y)P since a bias is intro- 
duced by its coupling to the classical particle coordinates. 
Some negative consequences of this fact are: i) The proof 
given in [ij that the total quantum-classical energy is con- 
served no longer applies. In fact energy is not conserved 
(no such claim was made in '^). ii) Unlike /, the distri- 



bution n(u, Ip; t) does not satisfy an evolution equation, 
that is, the knowledge of this distribution at some inter- 
mediate time t is insufficient to predict its future evolu- 
tion since the distribution of the Bohmian particle is also 
needed. Unfortunately, physics lies in n not /: If the 
proposal in the Letter is taken seriously as a prescription 
to evolve n from an initial time i to a final time t' (let 
us denote this operation by U{t' , t)) it will be found that 
[7(^2, ti)U{ti,to) differs from U{t2, to) since the proposed 
prescription requires to re-sample the Bohmian particle 
distribution at the intermediate time ti. 

A second problem is as follows. As is well known 
in pure quantum mechanics, given an ensemble of 
pure states, only the combination J2i Pi\''Pi) (''Pil 
evant to physics, since nothing else can be extracted 
from the knowledge of the expectation values of all 
possible observables. For the same reason, in the 
mixed quantum-classical case, only the quantity p{u) = 
J Dip n{u,ip) \'ilj)(tp\ (the mixed density matrix-density 
distribution in phase-space) characterizes the physical 
state. Note that there are many different n producing 
the same p. The trouble is that in the proposed method, 
two different initial distributions n corresponding to the 
same p will in general evolve into different p at later times 
(as it is readily shown), in other words, identical initial 
physical states p will branch into different states p{t) at 
later times, depending on which initial distribution n was 
chosen to represent it. This is a general problem when- 
ever the evolution is not directly based on the quantity 
p{u). 

This said, we have to distinguish between practical and 
theoretical points of view. At a practical level the method 
has been found to be useful and give good results in con- 
crete applications Our critique imphes, however, 
that at an exigent theoretical level the method fails to be 
fully internally consistent (this merely reflecting its in- 
trinsically approximate nature) and cannot be regarded 
as a final answer to the quantum-classical mixing prob- 
lem. 

This work is supported in part by funds provided by 
the Spanish DGICYT grant no. PB98-1367 and Junta 
de Andalucfa grant no. FQM0225. 



[1] O. V. Prezhdo and C. Brooksby, Pliys. Rev. Lett. 86, 3215 
(2001). 

[2] P. R. Holland, The quantum theory of motion (Cambridge 



2 



University Press., Cambridge, 1993). [4] E. Gindensperger, C. Meier, J. A. Beswick, and M.-C. 

[3] E. Gindensperger, G. Meier, and J. A. Beswick, J. Chem. Heitz, J. Ghem. Phys. 116, 10051 (2002). 

Phys. 113, 9369 (2000). 



